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METHOD ANI> SYSTEM FOR INDICATING ANISOTROPIC RESISTIVITY IN 

AN EARTH FORMATION 

1. BACKGROUND OF THE INVENTION 

1 . 1 Field of tlie Invention 

The present invention is generally directed to the analysis of underground earth 
formations, and, more particularly, to the determination of formation resistivity properties 
and/or profiles. 

1.2 Description of Related Art 

Knowledge of the resistivity properties of underground earth formations is 
fundamental in the evaluation and characterization of potential and existing hydrocarbon- 
bearing reservoirs. The terms "logging" and "log" are respectively used herein to 
designate any setjuential measurements and the recording of such measurements as a 
function of depth, and/or time of variations in a given characteristic of the formations 
around a borehok'. 

Formation resistivity anisotropy is a characteristic of reservoirs that can 
complicate formation evaluation. Many reservoir rocks exhibit resistivity anisotropy, 
especially when saturated with oil. There are several mechanisms, which can produce 
this anisotropy, among which are very thin sand-shale laminations, depositional changes 
in clean sandstons, and wind-distributed sands (aeolian formations). See Rubin, D. M., 
Cross bedding, bedforms, and paleocurrents. Society of Economic Paleontologists 
AND Mineralogists, Concepts in Sedimentology and Paleontology, 1 ; Klein et al., 
The petrophysics of electrically anisotropic reservoirs. Transactions of the SPWLA 
Thirty-Sixth Aiwual Logging Symposium, Paris, France, Jime 26-29, 1995, paper 
HH. Resistivity measurements from within a borehole drilled through these formations 
have been made; by a number of techniques employing the use of well-known 
instruments, including induction and laterolog (electrode) type devices. 
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In induction logging, a transmitter coil mounted on a well tool is lowered into a 
borehole and energized by an alternating current. Indications are then obtained of the 
influence of surrounding formations on the electromagnetic field established by the coil 
Usually such indications are obtained by observing the voltage induced in a receiver coil 
also mounted on the well tool in coaxial relationship v^th the transmitter coil and axially 
spaced apart therefrom. Conventional induction well tools include at least one transmitter 
and a receiver coil array. Such an induction tool is illustrated and described in U.S. Pat. 
No. 5,157,605, assigned to the present assignee. 

The laterolog technique comprises making a measurement with a well tool having 
a current emitting electrode and voltage electrodes one either side thereof which are 
operated to force current into the formation. The current is focused into the formations 
and returns locally, i.e., to the tool body. The formation resistivity is determined by 
measuring the cuirent flow. Conventional laterolog tools are illustrated and described in 
U.S. Pats. Nos. 5,754,050, 5,852,363, 6,023,168, and Documents EP-478409, EP- 
544583, all assigned to the present assignee. Document EP-544583 describes a laterolog 
tool including bcth a "deep" investigation mode and a "shallow" (LLS) investigation 
mode. 

The sensitivity of a resistivity logging tool to resistivity anisotropy depends on the 
physics of the measurement and the relative dip angle of the well bore of the formation. 
For example, 2MHz induction-type tools have no sensitivity to resistivity anisotropy 
when the formation layers are perpendicular to the well bore, but v^il be quite sensitive at 
relative dip angles above 60 degrees. See Luling et al. Processing and modeling 2 MHz 
resistivity tools in dipping, laminated, anisotropic formations^ Transactions of the 
SPWLA Thirty-Fifth Annual Logging Symposium, Tulsa, Oklahoma, June 19-22, 
1994, paper QQ. 

Induction logging tools are insensitive to resistivity anisotropy at low relative dip. 
The current density set up in the formation is limited to the horizontal direction. See 
Moran et al.. Effects of formation anisotropy on resistivity-logging measurements. 
Geophysics, V. 44, No. 7, July 1979, 1266-86; Anderson et al., The Response of 
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Induction Tools to Dipping Anisotropic Formations, Transactions of the SPWLA 
Thirty-Sixth Ai^jnual Logging Symposium, Paris, France, June 26-29, 1995, paper D. 
At higher relative dips the measured resistivity will combine both vertical and horizontal 
components, Hcwever, the effect does not depend on the coil array spacing, so the 
components are not resolvable. See Anderson et al, The Response of Induction Tools to 
Dipping Anisotropic Formations, Transactions of the SPWLA Thirty-Sixth Annual 
Logging Symposium, Paris, France, June 26-29. 

Laterolog logging tools have some sensitivity to the vertical component of 
resistivity. See Moran et al, Effects of formation anisotropy on resistivity-logging 
measurements, Gi^OPHYSics, V. 44, No. 7, July 1979, 1266-86; Chemali et al., The effect 
of shale anisotropy on focused resistivity devices. Transactions of the SPWLA Twenty- 
Eighth Annual Logging Symposium, London, England, June 30- July 2, 1987, paper H. 
However, with the dual-laterolog, the deep laterolog was not very sensitive to anisotropy. 
The shallow laterolog, which returns current to nearby electrodes, is more sensitive to 
anisotropy {See Moran et al, 1979), but invasion will reduce the effect of anisotropy. 

It is desii-able to obtain a reliable method and system for determining the 
resistivity properties and profiles of potential hydrocarbon-bearing zones in subsurface 
formations. Thus, there remains a need for improved methods and techniques to indicate 
formation resistivity anisotropy and evaluate the vertical component of resistivity. 

2. SUMMARY OF THE INVENTION 

Methods and systems are disclosed for combining induction tool and laterolog 
tool measurements to indicate that the formation resistivity is anisotropic and estimate the 
vertical component of resistivity even in the presence of invasion. 

One aspeci of the invention provides a method for determining a characteristic of 
a subsurface earth formation surrounding a borehole. The method comprises calculating 
at least one parameter representative of a property of the formation using log data 
obtained with a first well tool disposed within the borehole; modeling log data from the at 
least one calculated parameter; and comparing log data obtained with a second well tool 
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disposed within the borehole against the modeled log data to determine the formation 
characteristic. 

Another aspect of the formation provides a method for determining a 
characteristic of a subsurface earth formation surrounding a borehole. The method 
comprises disposing first and second well tools within the borehole to obtain 
measurement data; producing at least one parameter representative of a property of the 
formation from the measurement data obtained with the first well tool; producing 
theoretical data fi-om the at least one produced parameter; and comparing the theoretical 
data against the measurement data obtained with the second well tool to determine the 
formation characteristic. 

Another aspect of the invention provides a well logging system including first and 
second well tools adapted to be moveable through a borehole and an apparatus adapted to 
be coupled to the well tools, the apparatus being adapted to respond to data obtained with 
the well tools to determine a characteristic of a formation. The apparatus comprises 
means for calculating at least one parameter representative of a property of the formation 
using log data obtained with the first well tool disposed within the borehole; means for 
modeling log data from the at least one calculated parameter; and means for comparing 
log data obtained with the second well tool disposed within the borehole against the 
modeled log data to determine the formation characteristic. 

Yet another aspect of the invention provides a program storage device readable by 
a processor and eiicoding a program of instructions including instructions for performing 
operations in accord with the invention. 

3. BRIEF DESCRIPTION OF THE DRAWINGS 

Other aspects and advantages of the invention will become apparent upon reading 
the following detailed description and upon reference to the drawings in which: 

Figure 1 is a schematic diagram of the laterolog energization patterns obtained 
with a conventional laterolog tool. 
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Figure 2 is a schematic diagram of the current lines corresponding to the laterolog 
patterns of Figure 1. 

Figure 3 is a graph showing the response of the five laterolog arrays of Figures 1 
and 2 in a homogeneous anisotropic formation. 

Figures 4-7 are graphical illustrations of the responses of the five laterolog arrays 
of Figures 1 and 2 in invaded, anisotropic formations. 

Figure 8 is a graph showing the comparison of a laterolog tool response in an 
anisotropic formation vs. its response in an isotropic formation in accord with the 
invention. 

Figure 9 illustrates several plots of a ratio used in the calculations of the 
invention. 

Figure 10 illustrates the ratio plots of Figure 9 with an applied scaling factor. 

Figure 1 1 shows a crossplot for seventy ratios pertaining to cases of formation 
parameters in accord with the invention. 

Figure 12 shows an average value of a parameter for the ratios of Figure 11 in 
accord with the invention. 

Figure 13 is a graph showing a computed formation vertical resistivity component 
in accord with the invention. 

Figure 14 is a graph showing computed induction tool and laterolog tool logs in 
accord with the invention. 

Figure 15 is a graph showing the inversion of the induction tool data of Figure 14 
using a laterolog model in accord with the invention. 

Figure 16 is a graph showing the vertical resistivity component derived fi-om the 
data of Figures 14 and 15 in accord with the invention. 

Figure 17 shows field logs obtained from an induction tool and a laterolog tool in 
accord with the invention. 

Figure 1 8 shows the vertical resistivity component derived from the data of Figure 
17 in accord with the invention. 
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Figure 19 is a schematic diagram of a well logging system in accord with the 
invention. 

Figure 20 is a flow diagram of a method in accord with the invention. 
Figure 21 is a flow diagram of another method in accord with the invention. 

4. DEI AILED DESCRIPTION OF SPECIFIC EMBODIMENTS 

In the interest of clarity, not all features of actual implementation are described in 
this specification. It will be appreciated that although the development of any such actual 
implementation might be complex and time-consuming, it would nevertheless be a 
routine undertaking for those of ordinary skill in the art having the benefit of this 
disclosure. 

Figure 1 shows the energization pattems, referred to as "laterologs", obtainable 
from a laterolog-type tool 10 such as those disclosed in U.S. Pats. Nos. 5,754,050, 
5,852,363, 6,023,168, and Document EP-478409. The deepest array, mode 5, is 
considerably deepter than the LLS and should have considerable sensitivity to anisotropy. 
Figure 2 shows tlie current paths on several of the arrays of the tool 10. Note that the 
current flow lines turn parallel to the tool 10 axis and return locally, thereby increasing 
the sensitivity to the vertical resistivity component. 

Figure 3 shows the response of the tool 10 arrays (1-5) to an anisotropic 
homogeneous formation. The ratio of the apparent resistivity Rapp over the formation 
horizontal resisti\dty component Rth is plotted against the ratio of the formation vertical 
resistivity component Rty over Rth. This separation shows that the shallow arrays are 
much more sensii:ive to an anisotropic formation than the deeper arrays. Although the 
homogeneous medium anisotropy response looks favorable, what really matters is the 
response in the presence of invasion. 

Figures 4 through 7 show four representative cases of formation resistivity 
measurements obtained with a laterolog tool in an anisotropic formation having invasion 
out to a radius cf 60 inches. These figures show ratios of the formation horizontal 
resistivity component Rth over the formation invasion resistivity Rxo of 0.1, 0.2, 5, and 
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10. For all the crises Rtv ^ 10 x Rth. The invaded zone is isotropic with resistivity Rxo. 
In Figures 4 and 5, Rxo > Rth. In Figures 6 and 7, Rxo < Rth. For all these cases, the 
invasion response; confuses the anisotropy response. 

Figure 8 Sihows the anisotropic case of Figure 7 superimposed with the laterolog 
tool response to ;m isotropic virgin zone with formation resistivity Rt equal to Rth. As 
seen in Figure 8, the two curves representing Array 5 (RLA5) are almost parallel over the 
range of invasion radius out to 60 inches. The other shallower isotropic-vs.-anisotropic 
cases show separation at shallow invasion, but little or no separation at deeper invasion 
radii. The reasonable constancy of the difference between the isotropic and anisotropic 
Array 5 response suggests that the ratio of anisotropic to isotropic response or other ratios 
provide useful information. 

Designating Raa as the laterolog Array 5 log in an anisotropic formation, and Rai 
as the laterolog .\rray 5 log reading in an isotropic formation, the following ratio is 
derived: 



Figure 9 shows a plot of the above ratio for Rxo =1 and Rth =10 for X values from 1.4 to 
10, where 



It should be noted that the spread is almost proportional to X . This suggests that 



Rai 



{Raa - Rai) 




the ratio 



XRai 



Raa - Rai 



should be nearly constant. This ratio is shown in Figure 10. 
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If the ratio 



A = 



ARai 



(1) 



(Raa - Rai) 



is nearly constant for this case, then other values of Rxo and Rth should be looked at. 
Figure 1 1 shows a crossplot of 7 combinations of Rxo and Rth. There are 10 values of X 
ranging from 1.4 to 10. The spread in A is remarkably small. Figure 12 shovv^s the 
average value of A from the data of Figure 11. If an average value of A is taken to be a 
constant (or fitted function of the radius of invasion ri\ then the unknown X can be 
solved for 



One unknown remains in equation (2) - Rai, 

As mentioned above, induction tools have little or no sensitivity to Rtv in a 
vertical well. In such a case, the induction tool log data is inverted for Rxo, Rth, and ri 
using a lookup table, computer algorithm, model, or any other inversion means as known 
in the art. Also laiown in the art is the use of modeling techniques to simulate specific 
tool responses from given data. By applying a one-dimensional (ID) model of the 
laterolog tool response to the formation parameters obtained from the induction tool, 
Rxo, Rth, and r/, the laterolog log data in an isotropic formation {Rai) can be inferred. 

Then equation (2) can be solved for X , and then 



Figure 13 shows such a solution using a fitted value of ^ as a fimction of ru The values 
for A may be prec;omputed. The value for A ID model for the laterolog tool was derived 
using a simple a\^erage /-fimction, which is only a gross estimate for an example. The 
estimate of Rtv is within a factor of 2 of the actual value of Rtv in the uninvaded zone. 
Although a factor of 2 may be considered inaccurate for estimates of petrophysical 
quantities, such a(xuracy is ample for determining that anisotropy exists in a vertical well. 



X = A 



{Raa - Rai) 
Rai 



(2) 



Rtv= A^Rth. 



(3) 
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The accuracy may be improved by refining A, Both induction and laterolog tools have 
effects in 2D formations (layers and invasion in the layers) that are not accurately 
modeled v^ith ID' forward models. To test the usefulness of the ID models used in the 
recreation of Rai, it is advisable to try the algorithm in a 2D-formation model 

Figure 14 shows computed induction and laterolog tool logs in a simple 2D 
formation with 3 invaded, anisotropic layers of varying thickness. These layers have 
identical values of Rtv, Rth, Rxo, and rz. The thickness varies from 50 ft. to 10 ft. The 
induction tool log data reflects the receiver to transmitter spacings (10 in.- 90 in.) for the 
particular tool. 

Figure 15 shows a modeled laterolog tool log derived from the induction tool log 
data of Figure 14. The modeled log data represents the theoretical log values obtainable 
with the laterolog tool disposed within the isotropic formation. The modeled laterolog 
log data are then compared against the actual or apparent laterolog log data to compute 
Rtv as described above. Figure 16 shows the resulting Rtv log computed from the data 
of Figures 14 and 15. 

The Rtv log of Figure 16 shows anomalies in the non-invaded isotropic beds 
between the invaded anisotropic beds. As seen in Figure 14, the laterolog tool logs 
display shoulder effect in these beds and on either side of the first and third beds. Any 
source of error that leads to an increased separation between the laterolog tool logs 
predicted from the induction tool parameters and apparent laterolog tool logs will be 
interpreted as anisotropy. Overall, however, intelligent interpretation should remove 
most ambiguity, 

A field test was conducted comparing a laterolog tool against an induction tool in 
a formation with expected resistivity anisotropy. Figure 17 shows the induction tool and 
laterolog tool logs in a pay sand. Figure 18 shows the indicated vertical resistivity Rtv 
along with the induction tool Rth, Rxo, and ri inversion in a zone that exhibits the 
anisotropy signature. It also exhibits shoulder effect just as in the computed example of 
Figure 16. 
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4.1 Program Storage Device 

It will be apparent to those of ordinary skill having the benefit of this disclosure 
that the present invention may be implemented by programming one or more suitable 
general-purpose computers having appropriate hardware. The programming may be 
accomplished through the use of one or more program storage devices readable by the 
computer processor and encoding one or more programs of instructions executable by the 
computer for performing the operations described above. The program storage device 
may take the form of, e.g., one or more floppy disks; a CD ROM or other optical disk; a 
magnetic tape; a read-only memory chip (ROM); and other forms of the kind well-known 
in the art or subs€;quently developed. The program of instructions may be "object code," 
i.e., in binary form that is executable more-or-less directly by the computer; in "source 
code" that requires compilation or interpretation before execution; or in some 
intermediate form such as partially compiled code. The precise forms of the program 
storage device and of the encoding of instructions are immaterial here, 

4.2 Logging Slystem 

The present invention can be implemented in a logging system 200 with a first 
well tool 210 and a second well tool 215 adapted to be moveable through a borehole and 
an apparatus 220 coupled to the tools as shown in Figure 19. Certain conventional details 
are omitted in Figure 19 for clarity of illustration. The apparatus 220 comprises a surface 
computer (such as the general-purpose computer and program storage device described 
above) coupled tci the well tools 210, 215 by a wireline cable 225 or linked to the tool 
string's surface instrumentation as known in the art. 

The first v^ell tool 210 may be any laterolog-type tool that retums current locally, 
i.e., to the sonde body. The second well tool 215 may be any propagation or induction 
tool that provides multiple curves, preferably providing at least three independent 
measurements to solve for Rxo, Rt, and ri. In addition to the U.S. Patents mentioned 
above, other useal^le tools include the instruments described and illustrated in U.S. Pats. 
Nos. 5,339,037 and 4,968,940, both assigned to the present assignee. 
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The formation parameters can be determined in real-time by sending the 
measurement data to the surface as they are acquired, or it can be determined from a 
recorded-mode by recording the data on a suitable recordable medium. As known in the 
art, the measured data are transmitted from the tools 210, 215 to the surface computer 220 
by electronics (not shown) housed in the tools. The data may be sent to the surface 
computer 220 along the wireline cable 225 or by mud-pulse telemetry if implemented in a 
logging-while-drilling configuration. It will be understood that alternative means can be 
employed for communicating the acquired data to the surface as the precise form of 
communication is immaterial to the implementation of the disclosed techniques. 

Once received by the surface computer 220, the data can be recorded, processed, 
or computed as desired by the user to generate the formation resistivity profile. The 
profile can then be recorded on a suitable output record medium. Alternatively, some or 
all of the processing can be performed downhole and the data can be recorded uphole, 
downhole, or botti. It will be understood that the present invention can be implemented 
in conjunction with any suitable technique for keeping track of the tool 210, 215 depth 
within a borehole. 

Figure 20 illustrates a flow diagram of a method for determining a characteristic 
of a subsurface ibrmation surrounding a borehole in accord with the invention. The 
method comprises calculating at least one parameter representative of a property of the 
formation using log data obtained with a first well tool disposed within the borehole 300; 
modeling log data, from the at least one calculated parameter 305; and comparing log data 
obtained with a second well tool disposed within the borehole against the modeled log 
data to determine the formation characteristic 310. 

Figure 21 illustrates a flow diagram of a method for determining a characteristic 
of a subsurface formation surrounding a borehole in accord with the invention. The 
method comprises disposing first and second well tools within the borehole to obtain 
measurement datei 400; producing at least one parameter representative of a property of 
the formation from the measurement data obtained with the first well tool 405; producing 
theoretical data from the at least one produced parameter 410; and comparing the 
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produced theoretical data against the measurement data obtained with the second well 
tool to determine the formation characteristic 415. 

While the methods and apparatus of this invention have been described as specific 
embodiments, it will be apparent to those skilled in the art that variations may be applied 
to the structures and in the steps or in the sequence of steps of the methods described 
herein without departing from the concept and scope of the invention. All such similar 
variations apparent to those skilled in the art are deemed to be within this 
concept and sco])e of the invention as defined by the appended claims. 
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WHAT IS CLAIMED IS: 

1. A method for determining a characteristic of a subsurface earth formation 
surroimding a borehole, comprising: 

(a) calculating at least one parameter representative of a property of the 
foimation using log data obtained with a first well tool disposed within the 
bo rehole; 

(b) mctdeling log data from the at least one calculated parameter; and 

(c) comparing log data obtained with a second well tool disposed within the 
borehole against the modeled log data to determine the formation 
characteristic. 

2. The method of claim 1, wherein step (b) comprises modeling log data 
representative of log data theoretically obtainable with the second well tool 
disposed \^ithin the borehole. 

3. The method of claim 1, wherein step (c) comprises calculating a ratio between the 
log data ol^tained with the second well tool and the modeled log data. 

4. The method of claim 1, wherein step (c) comprises calculating the difference 
between the log data obtained with the second well tool and the modeled log data. 

5. The method of claim 4, wherein step (c) comprises multiplying the calculated 
difference by a predetermined factor. 

6. The method of claim 1 , wherein the first well tool is an induction-type tool. 

7. The methcd of claim 1 , wherein the second well tool is a laterolog-type tool. 

8. The method of claim 1, wherein the determined formation characteristic is a 
resistivity value or resistivity profile. 
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1 9. The method of claim 1, wherein the method is performed during or after drilling 

2 of the borehole. 

3 

4 10. A method for determining a characteristic of a subsurface earth formation 

5 surrounduig a borehole, comprising: 

6 (a) disposing first and second well tools within the borehole to obtain 

7 me;asurement data; 

8 (b) producing at least one parameter representative of a property of the 

9 formation from the measurement data obtained with the first well tool; 

10 (c) producing theoretical data from the at least one produced parameter of step 

11 (b); and 

12 (d) comparing the theoretical data of step (c) against the measurement data 

13 obtained with the second well tool to determine the formation 

14 chciracteristic. 

15 

16 11. The method of claim 10, wherein step (c) comprises producing theoretical data 

17 representative of measurement data obtainable with the second well tool disposed 

18 within the borehole. 

19 

20 12. The method of claim 10, wherein step (d) comprises calculating a ratio between 

21 the measurement data obtained with the second well tool and the produced 

22 theoretical data. 

23 

24 13. The method of claim 10, wherein step (d) comprises calculating the difference 

25 between the measurement data obtained with the second well tool and the 

26 produced theoretical data. 

27 

28 14. The method of claim 13, wherein step (d) comprises multiplying the calculated 

29 difference by a predetermined factor. 

30 
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1 15. The method of claim 10, wherein the first well tool is an induction-type tool. 

2 

3 16. The method of claim 1 0, wherein the second well tool is a laterolog-type tool 

4 

5 17. The method of claim 10 wherein the determined formation characteristic is a 

6 resistivity value or resistivity profile. 

7 

8 18. The method of claim 10, wherein the method is performed during or after drilling 

9 of the borehole. 

10 

11 19. A well logging system including first and second well tools adapted to be 

12 moveable through a borehole and an apparatus adapted to be coupled to the well 

13 tools, the apparatus being adapted to respond to data obtained with the well tools 

14 to determine a characteristic of a formation, the apparatus comprising: 

15 means for calculating at least one parameter representative of a property of the 

16 formation using log data obtained with the first well tool disposed within 

17 the borehole; 

18 means for modeling log data from the at least one calculated parameter; and 

19 means for comparing log data obtained with the second well tool disposed 

20 within the borehole against the modeled log data to determine the 

21 formation characteristic. 

22 

23 20. The system of claim 19, the modeling means further comprising means for 

24 modeling log data that is representative of log data theoretically obtainable with 

25 the second well tool disposed within the borehole. 

26 

27 21. The system of claim 19, the comparison means further comprising means for 

28 calculating a ratio between the log data obtained with the second well tool and the 

29 modeled log data. 
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1 22. The system of claim 19, the comparison means further comprising means for 

2 calculating the difference between the log data obtained with the second well tool 

3 and the modeled log data. 

4 

5 23, The system of claim 19, wherein the first well tool is an induction-type tool. 

6 

7 24. The system of claim 1 9, wherein the second well tool is a laterolog-type tool. 

8 

9 25. The system of claim 19, wherein the determined formation characteristic is a 

10 resistivity value or resistivity profile. 
11 

12 26. A program storage device readable by a processor and encoding a program of 

13 instructions including instructions for performing the apparatus operations recited 

14 in claim 19. 
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ABSTRACT 

A system and methods for indicating formation resistivity anisotropy and evaluating the 
vertical resistivity component by combining array induction measurements and array 
laterolog measurements. Formation measurements are obtained with an induction tool 
and a laterolog tool. Formation parameter data are derived from the induction tool 
measurements. The laterolog tool response is modeled using the derived formation 
parameter data. The formation resistivity anisotropy and/or the vertical resistivity 
component are determined from a comparison of the modeled laterolog tool response 
against the apparent measurements obtained with the laterolog tool. 
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